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Synchronization of many coupled oscillators is widely found in nature and has the potential to 
revolutionize timing technologies. Here we demonstrate synchronization in arrays of silicon nitride 
micromechanical oscillators coupled in an all-to-all configuration purely through an optical radiation 
field. We show that the phase noise of the synchronized oscillators can be improved by almost 10 
dB below the phase noise limit for each individual oscillator. These results open a practical route 
towards synchronized oscillator networks. 


Nano- and micromechanical oscillator arrays have the 
potential to enable high power and low noise integrated 
frequency sources that play a key role in sensing and the 
essential time keeping of modern technology [1-5 . The 
challenge with building scalable oscillator arrays is that 
micromechanical oscillators fabricated on a chip funda¬ 
mentally have a spread of mechanical frequencies due to 
unavoidable statistical variations in the fabrication pro¬ 
cess gim]. This dispersion in mechanical frequencies 
has a detrimental effect on the coherent operation in ar¬ 
rays of micromechanical oscillators. Here we show that 
arrays consisting of three, four and seven dissimilar mi¬ 
croscale optomechanical oscillators can be synchronized 
to oscillate in unison coupled purely through a common 
optical cavity field using less than a milliwatt of opti¬ 
cal power. We further demonstrate that the phase noise 
of the oscillation signal can be reduced by a factor of 
N below the thermomechanical phase noise limit of each 
individual oscillator as N oscillators are synchronized, 
in agreement with theoretical predictions 10, IT] . The 
highly efficient, low loss and controllable nature of light 
mediated coupling could put large scale nano- and mi¬ 
cromechanical oscillator networks in practice [T2HT8] . 

Synchronization is a ubiquitous phenomenon found in 
coupled oscillator systems mm- Heart beat is a result 
of synchronized motion of pace maker cells m , circadian 
rhythm arises because of coordinated body physiology 
[2Tj and global positioning system relies on synchronized 
operation of clocks. On the nanoscale, synchronization 
has been experimentally demonstrated in nanomechan¬ 
ical systems coupled through mechanical connections 
[3], electrical capacitors [9], off-chip connections [6 and 
through an optical cavity mm- However, these demon¬ 
strations were limited to only two oscillators. Achieving 
synchronization in large micromechanical oscillator net¬ 
works requires scalable oscillator units and efficient and 
controllable coupling mechanisms mm in]. 

Here we experimentally demonstrate that arrays of 
free running micromechanical oscillators can be synchro¬ 
nized when coupled purely through a common electro- 



FIG. 1: Concept and devices (a) Concept of mediating cou¬ 
pling between mechanical oscillators (yellow) through a global 
optical field (blue). The optical field provides energy for each 
mechanical oscillator to vibrate at their natural frequencies 
Qij and also provide coupling between each mechanical oscil¬ 
lator forming an all-to-all coupling topology. When the opti¬ 
cal coupling is strong, the oscillators synchronize and vibrate 
at a common frequency, (b) A schematic of each individ¬ 
ual double-disk. The edges are partly suspended to allow for 
mechanical vibration (c) Cross section of a double-disk show¬ 
ing the mechanical and the optical mode shapes, (d) Optical 
microscope images of coupled optomechanical double-disk os¬ 
cillator arrays. The oscillators are mechanically separated by 
a narrow gap (~ 150 nm) and coupled solely through the op¬ 
tical evanescent field. The squares and strings are support 
structures for tapered optical fibers. 


magnetic field as predicted by theories [12] 13 . A con¬ 
ceptual view of an array of mechanical resonators cou¬ 
pled through light is illustrated in Figure [l^i. Each op¬ 
tomechanical oscillator (OMO) possesses a slightly differ¬ 
ent frequency of mechanical oscillation (f^) and is only 
connected through a common optical field (blue back¬ 
ground). When a continuous wave laser is coupled to 
a common electromagnetic field mode spanning several 
micromechanical oscillators, the light can provide both 
the drive for self-sustaining oscillations and the necessary 
coupling between the individual oscillators for synchro¬ 
nization through optical forces. When the laser power 
is just above the self-sustaining oscillation threshold of 
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the mechanical oscillators, they are expected to vibrate 
at their natural frequencies fV When the laser power 
is high so that the optically mediated coupling is strong 
enough to overcome the difference in the mechanical 
oscillators can reach synchronization. 

The effective coupling between the mechanical res¬ 
onators can be visualized through the following equation 

Xi + TiXi + H? Xi = F c ( pt 

dpt a \b(xi,...,Xj )\ 2 

where XiXifii are the mechanical displacement, damp¬ 
ing and the mechanical frequency of the ith OMO and 
is the amplitude of the coupled optical su¬ 
permode that spatially spans all cavities in the array. It 
is clear from the equation above that the optical force 
(F opt ) depends on the energy stored in the optical su¬ 
permode which is affected by the displacement of each 
individual cavity. Therefore the optical field provides an 
effective nonlinear mechanical coupling between the dif¬ 
ferent oscillators that form the basis for synchronization 
in m ns]. The onset of synchronization, which intrinsi¬ 
cally relies on nonlinearity [23], could therefore be cap¬ 
tured as (Topt) is increased through increasing the optical 
driving power [7]. 

The individual oscillator we use is a double-disk OMO 
(Figures [lj) and [1]:) composed of two free-standing silicon 
nitride circular edges that support high quality ( Q ) factor 
optical and mechanical modes Om]. The co-localized 
modes shown in Figure [TJ^ lead to a strong coupling be¬ 
tween the optical and the mechanical degree of freedom. 
When the cavity is excited by a continuous wave laser 
above the oscillation threshold, the free-standing edges 
oscillate coherently and modulate the laser producing a 
radio frequency (RF) tone at the mechanical frequency 
of the vibrating edges. Fabrication variation causes the 
mechanical frequency of different OMOs in our arrays to 
spread around ±1 MHz centered at 132.5 MHz [26] . 

We fabricate micromechanical oscillator arrays with 
double-disk OMOs that are optically coupled through the 
evanescent field. The OMOs are physically separated by 
a narrow gap (^150 nm) which precludes any mechani¬ 
cal connections while the optical evanescent field can still 
propagate through the gap. Mechanical coupling through 
the substrate connection is negligible as the mechanical 
mode we excite is a high Q mode that is well isolated 
from the substrate [24]. 

We excite the optical supermode that spatially span 
over all cavities to ensure that there is optical coupling 
among all cavities (Fig. El dashed boxes). The strong 
optical coupling between the optical modes of each in¬ 
dividual cavity leads to the formation of optical su¬ 
permodes b m that have different optical frequencies and 
spatial geometries [26] . Figure El illustrates the spatial 
intensity profile of different optical supermodes b m when 
the optical resonant frequency of individual cavity uoi is 



b 



FIG. 2: Experimental configuration (a) Optical super¬ 
modes spatial structures. The colored halos show where the 
optical cavity field resides for different types of arrays. The 
more opaque colors illustrate higher cavity field intensities 
when compared to the rest of the cavities. The supermodes 
that spatially span over all cavities with equal intensities are 
identified by dashed lines, (b) Experimental setup. The 
coupled optomechanical oscillator array is placed in a vac¬ 
uum chamber and excited by a tunable infra-red (IR) cam¬ 
era through a tapered optical fiber. The optical power and 
polarization are controlled by an variable optical attenuator 
(VOA) and a fiber polarization controller (PC). The optical 
transmission is detected by an amplified photodiode (PD) and 
analyzed by an oscilloscope and a spectrum analyzer. 


identical. The higher intensity regions are illustrated by 
higher opacity of the halos around the cavities. We posi¬ 
tion a tapered optical fiber to the close proximity of one 
OMO in the arrays to couple light to the spatial evenly 
distributed optical modes (dashed lines in Fig. [2^i) while 
using an infrared (IR)-camera to monitor the scattered 
intensity from the arrays making sure all OMOs are ex¬ 
cited. We monitor the transmission through the tapered 
fiber by an amplified photodiode and feed the electrical 
signal to a spectrum analyzer. 

We show the onset of synchronization by increasing 
the excitation laser power which effectively increases the 
coupling between the OMOs. The laser wavelength is 
blue detuned relative to the resonance of the optical 
supermode that evenly spans all the OMOs (Fig. El 
dashed boxes) enabling optomechanical amplification. In 
the three coupled OMO array, as the laser power in¬ 
creases well beyond the oscillation threshold for each in¬ 
dividual oscillator, the radio frequency (RF) spectrum 
of the OMOs show many strong oscillation peaks and 
a broad noise floor (Fig. |3^i) . The distinct oscillation 
peaks form because Qi is different for each OMO and 
they beat to generate many RF tones [271128]. The in¬ 
crease in the noise floor is likely due to finite interaction 
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FIG. 3: Synchronization in arrays of OMOs Optical power spectrum of three- (a) four- (b) and seven- (c) OMO system as 
the input optical power increases. The vertical scale is from -110 dBm to 0 dBm for each trace. Synchronization is characterized 
by the sudden noise floor drop and the emergence of a single frequency in the optical power spectrum as indicated in the graphs. 
The disorder in natural mechanical frequencies and incoherent dynamics before the onset of synchronization is evident from the 
RF peaks and the broad noise floor. The seven resonator array (c) shows multiple changes of noise shapes before eventually 
synchronize indicating the presence of multiple oscillation states as a result of many OMOs. 


between the mechanical modes mediated by the optical 
field but not yet strong enough to transition into a locked 
state 01221 EDI. As the laser power further increases to 
Pin = 280 /iW, the onset of synchronization (Fig. Hf) is 
evident as the peaks on the RF spectrum merge into a 
single large peak and the noise floor is reduced. The much 
weaker sidebands around the main oscillation signal are 
due to the much weaker oscillatory motion induced by 
thermal force displacing the OMOs from the synchro¬ 
nized state 0. In four and seven coupled OMO arrays, 
similar to the three-cavity system, we observe beating 
between different mechanical modes and a broad noise 
floor when the optical power is below the synchroniza¬ 
tion threshold. As the laser power is increased, a single 
oscillation peak appears accompanied with a sudden drop 
in the noise floor signifying the onset of synchronization 
(Fig. |pc). 

We show that in large arrays of OMOs, the phase 
noise of the synchronized signal can be reduced below 
the thermomechanical noise limit of an individual OMO 
by almost 10 dB. The phase noise of the modulated out¬ 
put light is expected to drop as the oscillators are syn¬ 
chronized 0 ng mj. We measure the phase noise of 
our oscillators at 10 kHz offset from the carrier oscilla¬ 
tion frequency, where the phase noise of our oscillator is 
dominated by thermomechanical fluctuation [26. 13HI331 . 
a fundamental limit imposed to the mechanical oscillator 
due to the thermal bath of the environment. In figure 
HP we show the measured phase noise in a 1 x 2 OMO 
array [8]. As shown in figure [4^, the single OMO phase 


noise [26] at low optical power is ^ —60 dBc/Hz and 
gradually increases as the laser power is increased. The 
increase of phase noise is due to phase slipping between 
the two OMOs [6]. As the coupling between the OMOs 
increases with increasing laser power, they synchronize. 
As expected, we observe the phase noise drops by ~ 3 
dB as the two OMOs move from one OMO oscillating 
state to a synchronized oscillation state. Since the oscil¬ 
lators are nearly identical, synchronized oscillations can 
be viewed as two oscillators operating coherently provid¬ 
ing a larger effective mass while not reducing the oscilla¬ 
tion frequency [32]. In figure [4^, we show the measured 
phase noise of each large array of oscillators by driving 
the system at high optical powers at the optimal opti¬ 
cal detuning where the phase noise is a minimum [26] . 
The lowest phase noise measured in each array of dif¬ 
ferent sizes is plotted in figure [4 ]d. The measured phase 
noise follows the 1/N dependence predicted by theory 

mmmm- 

The drop in phase noise can also be used to determine 
the number of synchronized OMOs in a single array os¬ 
cillating in different states. We measure the phase noise 
in the 2x2 array as the oscillators change from a state 
where only two OMOs are oscillating to a state where 
all four OMOs are oscillating, as we infer from the light 
scattering intensities captured on the IR camera. Figure 
H shows the power spectrum of the transmitted light 
when the laser is tuned from an optical mode that spans 
two cavities to an optical mode that spans all four cav¬ 
ities (Fig. IT) while staying at the same optical power. 
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FIG. 4: Phase noise in synchronized arrays (a) phase noise in a two-OMO system at 10 kHz carrier offset as a the laser 
power is increased. The noise increases due to mode competition between possible oscillation states and then decreases by ~ 3 
dB below the noise level of one OMO oscillation state, (b) The phase noise of the synchronized oscillation signal for different 
sizes of OMO arrays. The grey curve is the phase noise level predicted by theory for near-identical synchronized oscillators (c) 
Power spectrum of a state where four OMO are oscillating (black) and of a state two OMO oscillating (purple). The phase 
noise drops by ~ 3 dB following the transition. The inset shows a zoom in near the peak and captured IR images for the two 
and four OMO synchronized state respectively. 


Following the transition, the four OMO oscillation state 
shows an increase of ^ 3 dB in the oscillation signal and 
^ 3 dB drop in the phase noise. At the same time, all 
four resonators light up on the IR camera. The drop in 
phase noise and the change of scattering intensity on the 
IR camera strongly indicate that the array changes from 
two to four synchronized oscillators. 

In conclusion, we demonstrate synchronization in in¬ 
tegrated arrays of micromechanical oscillators coupled 
through a common optical field. We show the onset of 
synchronization in the arrays by tracking the emergence 
of a single oscillation frequency in the optical power spec¬ 
trum. Synchronization is further corroborated by the 
drop in phase noise in the oscillator arrays. The reduc¬ 
tion of phase noise with oscillator array size and the seal- 
ability of our devices could enable low noise and high 
power integrated frequency sources. Our work paves a 
path towards large scale monolithically fabricated oscil¬ 
lator networks that have the potential to compete with 
the performance of bulk resonators and to exhibit rich 
nonlinear dynamics opening the door to novel metrology, 
communication and computing techniques [35 , 36]. 
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SUPPLEMENTAL INFORMATION 
Optical Transmission 



FIG. SI: Mode excitation geometry (a) In a square lattice, only one direction of the optical mode is excited, (b) In a triangular 
lattice, the coupling geometry determines that as the light from the waveguide excites the clockwise mode of OMO 1 which in 
turn excites the counter-clockwise mode of OMO 2 and OMO 3. However, the counter-clockwise mode of OMO 3 also excite 
the clockwise mode of OMO 2 (blue arrow), which therefore excites all the counter propagating wave of the entire system. 


The challenge in optically coupling large arrays of mechanical structures is the variation of the optical frequencies 
of each individual OMO due to variation in the exact fabricated dimensions. This can be in principle compensated by 
integrating optical tuning mechanisms such as heaters. In our device, we overcome this challenge by ensuring that the 
optical coupling is large compare to the frequency variation and the individual decay rate of each cavity. The optical 
coupling strength we designed (ftfaj) ~ 5 GHz) is much larger than the optical decay rate of the cavity ( 7 \ ~ 400 
MHz) and to the cavity optical frequency spread from fabrication variations (A uj < 2 GHz). This ensures that the 
optical modes are strongly coupled which means the light travels between the cavities many times before it is lost via 
other channels. 

We estimate the spatial structure of the optical supermodes using coupled mode theory m In the 2 x 2 OMO 
array, the coupled modes can be solved essentially using the same method as described in previous work of two coupled 
double disk OMOs [ 8 ]. The evolution and coupling of each individual optical mode is described by a matrix where 
its diagonal entries describe the eigenfrequencies and the off-diagonal entries describe the coupling between these 
modes. In the case for the three- and seven- OMO arrays arranged in a triangular lattice, we must consider both the 
clock-wise and counter-clock wise optical modes, as both will be excited through a laser coupled to one direction even 
without any backscatters (Fig. SI). 

The coupled mode equation that governs the optical part reads 
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where 7 i,u)i are the total damping rate and the resonant frequency of each individual optical mode is the 

optical coupling strength between ith and jth OMO. \s(t )\ 2 is the optical power that excites the array of resonators 
through the first cavity. r) c is the optical criticality factor, the ratio between the cavity damping rate and the external 
coupling rate. For triangular lattice cases, the matrix dimension is increased from n to 2 n where the both the cw 
and ccw degrees of freedoms are taken into account. In the strong coupling regime (/^j >> 7 i individual 

optical modes hybridize and fully splits into optical supermodes b m where some supermodes span over all cavities. 
The spatial geometry of these supermodes can be analyzed through rewriting equation [Si] in the optical supermode 

N 

basis b m = c i a i where <7 indicates the relative optical field amplitude in cavity i. To conveniently visualize the 

i =1 
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FIG. S2: The optical transmission in the (a) three, (b)four and (c) seven cavity arrays respectively. The well split resonant 
mode and the high extinction of each split supermodes shows that the optical modes are strongly coupled. 



FIG. S3: Image of the scattered light when sychronization is onset for (a) three-cavity (b) four-cavity and (c) seven-cavity 
arrays. The scale bars are 20 /xm. 


spatial geometry of the optical supermode, we approximate 7 * and f \ of individual OMO being identical. We can 
then exactly diagonalize the coupling matrix in equation [Si] through Jordan matrix decomposition to find the new 
optical supermodes basis in superposition of the individual optical mode a$. The geometric representation of the 
superposition is shown in figure 2 in the main text. 

Figure [S3] shows images from the IR camera when the supermodes spanning over all cavities are excited. It is clear 
from the image that all cavities light up with approximately equal intensities. The strong coupling is also visualized 
by the fully split optical supermodes shown in the optical transmission in figure |S 2 | 


Mechanical modes 


The mechanical modes of the double disk are the flapping mechanical mode where we can control the mechanical 
frequency of the mode by the undercut depth of the sacrificial silicon dioxide layer. In the particular device we used in 
the experiment, we fabricated the devices with mechanical frequencies near 132.5 MHz and quality factor near 1000 
(Fig. S4). Typical optical quality factor of these devices is ~ 500,000. interaction strength between the optical and 
the mechanical mode is characterized by the optomechanical coupling strength g om = ^ = 2tt x 50 GHz/nm for our 
devices 0. 


In the OMO arrays, when the system is coupled to a laser at very low power, we measure a distribution of mechanical 
frequencies corresponding to different OMOs. Figure |S5| shows the mechanical spectrum in the three, four and seven 
cavity arrays respectively. In the three OMO array, the mechanical frequencies are clearly different whereas in the 
four and seven OMO arrays, the mechanical mode frequencies are different but also overlapping. 
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FIG. S4: (a) Thermal mechanical spectrum of a single doubled-disk 0M0 forming the array in this experiment, (b) Self- 
sustaining oscillation spectrum of a single double-disk OMO. 





Mechanical Frequency (MHz) Mechanical Frequency (MHz) Mechanical Frequency (MHz) 

FIG. S5: The mechanical spectrum in the (a) three, (b)four and (c) seven cavity arrays at low optical powers respectively. The 
disorder in mechanical frequencies is evident from the splitting of the peaks. 


Phase Noise 

We measure the phase noise of our oscillator by feeding the output of the photodiode (New Focus 1811) to a 
spectrum analyzer (Agilent E4407B). To measure the phase noise at 10 kHz offset frequency, we utilize the spot 
frequency measurement option in the phase noise measurement module (Option 226) of the spectrum analyzer which 
enables fast noise spectral power measurement at a single offset frequency. 

In synchronized oscillator arrays, the phase noise reduces as the number of synchronized oscillators increase. Ex¬ 
periments and theories mm show that the noise reduces as 1 /N where N is the number of synchronized oscillators. 
This can be intuitively understood as coordinated oscillators while having the same frequency but having an increased 
effective mass which improves their resistance to thermo mechanical fluctuations. In our system, proving synchroniza¬ 
tion by studying the dynamics of individual cavities in the arrays is desirable but it is technically challenging given 
the compact size of the array. Therefore we choose the phase noise as a figure of merit as it is fundamentally limited 
by the thermomechanical process. 

The phase noise performance of an optomechanical oscillator, when limited by thermomechanical noise, improves 
when the laser power is increased. At high laser driving powers, the phase noise roughly stays constant as the laser 
power is further increased. We calibrate the single oscillator phase noise performance by measuring the phase noise 
at 10 kHz carrier offset at different laser driving powers. Figure [S6] clear shows that the phase noise does not vary 
significantly when the power is above 4 x P t h • where P t h = 10 /iW is the typical threshold power for a single OMO. 
This is in agreement with the experimental observation reported by Mani et. al [33]. and theory by Fong et. al m- 

Leesons model [38] predicts the phase noise of self-sustaining oscillators from the linewidth of the oscillation signal. 
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FIG. S6: Spot frequency phase noise of a single OMO at 10 kHz offset from the carrier frequency as a function of laser power. 
Pth = 10 fiW is the threshold power for a single OMO to oscillate. 



FIG. S7: FEM simulation of the mechanical mode, (a) Geometry of the structure, (b) Simulated mechanical displacement 
profile. The antisymmetric mechanical mode simulated has a vibration frequency of 131 MHz with an effective mass of 70 pg. 

The linewidth (AO) of a thermally limited oscillator is given by 

Aft = (S2) 

where ks is Boltzmann’s constant, T is the temperature, P out is the output power of the oscillator and T m is the 
natural damping rate of the oscillator. Since the double disk we use have a small effective mass m e ff (Fig. |S7| ), which 
means a low oscillator power, the phase noise is limited by thermomechanical noise. The oscillator power P out is given 
by, 


PDut — 2 meff ^ ^ F r 


(S3) 


where x is the maximum displacement amplitude of the oscillator. Equation (S2) can be represented as 


AO _ ^ b T _ 

m eff Q 2 x 2 71 


(S4) 


In our experiment, since the offset frequency A/ we are interested in is much larger than the linewidth (A Q) of our 
oscillator, the phase noise L(Af) in dBc/Hz predicted by Leeson’s eauation|38l can be simplified as, 
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HAf) = 10Log„ (SB) 

Substituting the parameters (T = 300 K; m e q = 70 pg; T m = 2tt x 110 kHz) and x ~ g om /7 = 10 pm to equation 
we obtain a theoretical estimate of the thermo mechanical noise limit of our oscillators which is —60 dBc/Hz at 10 kHz 
carrier frequency offset in good agreement with our experimental measurements and previous frequency dependent 
phase noise measurements [M. 

We expect the laser noise to be a minor contribution to the phase noise of our OMOs at 10 kHz carrier offset. 
Although diode lasers typically have relatively high phase noise, at this frequency, the phase noise of the laser is not 
relevant because the optical force is only affected by laser intensity. The cavity decay rate ^ 400 MHz) is much 
higher than 10 kHz which means the laser phase noise does not get converted in to intensity noise due to the cavity 
response. Typical relative intensity noise (RIN) of external cavity tunable laser at low frequencies is below -100 
dBc/Hz [315]. We measure the RIN of the laser we use (Tunics Reference) at 50 /iW to be -112 dBc/Hz at 10 kHz 
offset through the photodiode (New Focus 1811) and the spectrum analyzer (Agilent E4407B). As laser shot noise is 
much lower than technical noise at this laser power and frequency, we expect minimal contribution from the shot 
noise of the laser. Our analysis is consistent with previous theory m- 
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